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ABSTRACT 

Direct detection of extrasolar planets will require imaging systems capable of un- 
precedented contrast. Apodized pupils provide an attractive way to achieve such con- 
trast but they are difficult, perhaps impossible, to manufacture to the required toler- 
ance and they absorb about 90% of the light in order to create the apodization, which 
of course lengthens the exposure times needed for planet detection. A recently pro- 
posed alternative is to use two mirrors to accomplish the apodization. With such a 
system, no light is lost. In this paper, we provide a careful mathematical analysis, 
using one dimensional mirrors, of the on-axis and off-axis performance of such a two- 
mirror apodization system. There appear to be advantages and disadvantages to this 
approach. In addition to not losing any light, we show that the nonuniformity of the 
apodization implies an extra magnification of off-axis sources and thereby makes it 
possible to build a real system with about half the aperture that one would otherwise 
require or, equivalently, resolve planets at about half the angular separation as one can 
achieve with standard apodization. More specifically, ignoring pointing error and stel- 
lar disk size, a planet at 1.7 X/D ought to be at the edge of detectability. However, 
we show that the non-zero size of a stellar disk pushes the threshold for high-contrast 
so that a planet must be at least 2.5X/D from its star to be detectable. The off-axis 
analysis of two-dimensional mirrors is left for future study. 



Subject headings: Extrasolar planets, coronagraphy, point spread function, apodiza- 
tion 
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1. Introduction 

In a recent paper, Goncharov et al. (2002) propose using two mirrors to convert a uniform- 
intensity beam to one of nonuniform intensity — that is, an apodized exit beam. The motivation 
for their work was to "mask" the secondary mirror in designs of large two-mirror telescopes with 
fast spherical primaries. As an application of the idea, they point to radio and submillimeter- 
wave telescopes since for these it is important to suppress sidelobes and reduce background noise. 
However, its applicability to high-contrast imaging in the context of the Terrestrial Planet Finder 
(TPF) project immediately caught our attention and the result is the work described in this paper. 

Independently of Goncharov et. al., Guyon (2003) studied the same two-mirror apodiza- 
tion idea and analyzed its applicability to extra- solar planet finding. His work was motivated by 
Labeyrie's paper on pupil densification (Labeyrie (1996)). The idea behind pupil densification 
is to reduce the intensity of the diffraction side lobes by remapping the sparse input pupil of an 
interferometer to a more compact exit pupil. Pupil densification provides the advantage of a better- 
concentrated image at the cost of field of view, since it violates the "golden rule" of combining 
beams (Traub (1986)); this rule states that, to provide a large field of view, the exit pupil of an 
optical system must be a scaled version of the input pupil. Since pupil densification is a special 
case of pupil remapping, both will suffer field of view loss if the scaling factor is not constant. 

Nevertheless, the potential gains from pupil mapping are so significant that it is worthwhile 
exploring this technique for the purpose of extrasolar planet searching, and this is the goal of the 
present paper. 

In his paper, Guyon argues that the two-mirror approach provides a system with (i) a smaller 
inner working angle (IX/ D vs. AX/ D) and consequently a smaller overall aperture and (ii) a 
higher throughput with a corresponding reduction in exposure time. Here, and throughout the 
paper, A is wavelength and D is the diameter of the primary mirror. Guyon also addresses, at least 
qualitatively, the issue of off-axis performance of the system. 

In this paper, we give a careful mathematical analysis of the on-axis and off-axis performance 
of this two-mirror system. We deal exclusively with one-dimensional mirrors in two-dimensional 
space. While the on-axis analysis is easily extended to real two-dimensional mirrors in three space 
dimensions, the off-axis analysis, which is the crux of this paper, is much more difficult in this 
higher dimension. Hence, we leave such analysis to future work. 

The first and most obvious advantage of the two-mirror system is that throughput is increased 
since no light is lost. Furthermore, an interesting and unexpected advantage of this system is that 
there is an inherent magnification introduced by a nonuniform apodization. For example, a planet 
at 5X/D in the sky actually appears at about 12X/D in the image plane. This improves the effective 
inner working angle by a factor of 2.4 and explains, we believe, Guyon's claims about improved 
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inner working angle. Consequently, such a system for planet finding could be implemented with 
a smaller mirror than would otherwise be required or, equivalently, could detect planets less than 
half as far from their parent star with the same aperture. 

There is, however, a downside. We show that the inner working angle of an ideal system 
(with an on-axis source) is still AX/ D and furthermore that even tiny off-axis light sources, such as 
arising from the nonzero size of the stellar disk, increase the inner working angle to at least 6A / D. 
So, a star at 2.5X/D in the sky will appear at 6X/D in the image plane and hence lie on the edge of 
detectability. This still represents a significant improvement over the AX/D, which is the best one 
can do with standard apodization or with pupil masks. 

In the next section, we derive the basic equations relating the shapes of two one-dimensional 
mirrors to the apodization they produce. Then, in Section 3, we give a careful analysis of the off- 
axis performance of this two-mirror system. Section 4 then goes through some specific examples, 
including the case of two parabolas, where everything can be computed explicitly. Finally, in 
Section 5, we show the results of numerical computations using the known optimal apodizer. 

We end the present section with a word about our notation. Given any function of a single 
real variable, we use a prime to denote the derivative of this function with respect to its variable. 
So, for example, h'(x) = dh(x)/dx. 



2. Mirror-Based Apodization 

Figure 1 illustrates the use of two mirrors to produce an apodization. Parallel light rays come 
down from above, reflect off the bottom mirror (on the left), bounce up to the top mirror (on the 
right), and then exit downward as a parallel bundle with a concentration of rays in the center of the 
bundle and thinning out toward the edges — that is, the exit bundle is apodized. A focusing element 
on the x-axis then creates an image. 

Let h(x), a < x < b, denote the height-function defining the lower mirror and let g(x), 
— f < x < f , denote the height-function defining the upper mirror. We suppose there is a one-to- 
one correspondence between x's in the interval [a, b] and the interval [— f , f] that is determined by 
the desired apodization. Let x(x) denote the value in [— f , f] associated with x in [a, b}. We assume 
that the function x, which we call the transfer function, is strictly monotone (either increasing or 
decreasing). It then has a functional inverse. Let x(x) denote this inverse function. The apodization 
A(-) associated with this correspondence is given by 

— (x) = A(x). (1) 
Adding a boundary condition, such as x(— f ) = a, makes it possible to derive the correspondence 
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from the apodization. 

The shapes of the two mirrors are determined by enforcing equality between angle-of-incidence 
and angle-of-refiection (i.e., Snell's law) at the two mirror surfaces. The tangent vector to the 
first mirror is (l,h'(x)). The unit incidence vector is (0,1) and the unit reflection vector is 

(x — x, g(x) — h(x))/S(x, x), where 



S(x, x) = yjjjb — x) 2 + (g(x) — h(x)) 2 (2) 

denotes the distance between the point (x, h(x)) on the first mirror and (x,g(x)) on the second 
mirror. The scalar product between the tangent vector and the unit incidence vector must be equal 
to the negative of the scalar product between the tangent vector and the unit reflection vector. 
Equating these inner products, we get the following differential equation for h: 

x — X 

= 77, rr 7-77 7^7- (3) 

S(x,x) - h(x) + g(x) 

Applying the same technique to the surface of the second mirror at x, we get the same expression 
for g' and so 

g'(x) = h'(x). (4) 

Equation (4) says that, for a given ray, the slopes of the mirrors are locally the same, which is to 
be expected since the input ray is parallel to the output ray. From equations (3) and (4), it appears 
that the differential equations for h and g are a coupled system of ordinary differential equations. 
Our first task is to show how to decouple them. 

Let 

P (x, x) = S(x, x) +g(x) - h(x), (5) 

denote the denominator in (3). It is the total optical path length from a point on any given hori- 
zontal line drawn above the first mirror (and below the second) vertically down to the first mirror 
diagonally up to the second mirror and finally vertically back down to the same horizontal line 
below the second mirror. 



Theorem 1 The optical path length P is a conserved quantity for on-axis rays. 



Proof. If we regard x as a function of x and differentiate with respect to x, we get 

^ (S(x, x) + g(x) - h(x)) = G(x) - H(x)^, (6) 
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where 

g(x) = 

H(x) = ^- X + ( ^g- fe(x))fe/(:E) +^). (8) 

From (4), we see that G(x) = H(x). Using (3) to eliminate h'(x) from (8), we get after some 
straightforward algebraic manipulation that H(x) = and therefore P is a constant. □ 

Henceforth, we write P for this optical-path-length invariant. 

Corollary 1 The differential equations for h and g are decoupled: 

h'(x) = g'(x) = (9) 

-TO 

Isolating S(x, x) on one side of the equation in (5) and squaring, we can solve for the difference 
between g(x) and h(x): 

m - Kx) = + f . do) 

Hence, it is only necessary to solve one of the differential equations. The other mirror surface can 
then be computed from this difference. 



3. Off -Axis Performance 



Suppose now that coherent light enters at angle 6 from on-axis as shown in Figure 2. The light 
that strikes at position (x, h(x)) on the first mirror reflects off at a different angle and strikes the 
second mirror at a new position shifted by Ax from the on-axis strike point: (x + Arc, g(x + Ax)). 
Light then reflects off the second mirror at angle 9 and intersects the rr-axis at position x. We need 
to compute how Ax, 9, and x depend on 9 with the ultimate goal of using these expansions to 
compute the image-plane point-spread function as a function of 9. 



Lemma 1 Ax = S(x, x) 9 + o{9). 
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Proof. The unit incidence vector is (— sin 9, cos 9) and the unit reflection vector is (x + Ax — 
x, g(x + Ax) — h(x))/S(x, x + Ax). As before, Snell's law can be expressed by equating the inner 
product between the unit incidence vector and the mirror's tangent vector with the negative of the 
corresponding inner product for the unit reflection vector. We get that 

• x a (x + Ax - x + h'(x)(g(x + Ax) - h(x))) 
— smO + h (x) cosU = — — . (11) 

v ' S(x,x + Ax) 

Linearizing the left-hand side in 9 and the right-hand side in Ax and noting that the constant terms 
in the two linearizations match, since they give the on-axis equation satisfied by h'(x), we get 

(\ + h'(x)q'(x) x — x + h'(x)(q(x) — h(x)) . ,, ...\ „ ns 

9 ~ of V q( 3 — (x-x + g'(x)(g(x) - h(x))) Ax. (12) 

\ b{x, x) b{x,x)- i J 

Using Corollary 1, the first term on the right simplifies to (1 + h'(x) 2 )/S(x,x) and the second 
term simplifies to h'(x) 2 /S(x, x). Hence, we get that to first order 9 = Ax/S(x, x) from which 
the lemma follows. □ 



Lemma 2 9= Ax + o(Ax) = A(x)9 + o(9). 

S(x, x) 



Proof. The proof is similar to that of the previous lemma. The unit reflection vector is (sin 9, — cos 9) 
and the unit incidence vector is just the negative of what it was before: — (x + Ax — x, g(x + Ax) — 
h(x))/ S(x, x + Ax). The mirror's tangent vector at x + Ax is (1, g'(x + Ax)). As before, Snell's 
law can be expressed by equating the inner product between the unit incidence vector and the 
mirror's tangent vector with the negative of the corresponding inner product for the unit reflection 
vector. We get that 

. ~ ^ (x + Ax - x + g'(x + Ax)(g(x + Ax) - h(x))) . ^ _ > 

*»»-!f(*+**)<*»9 = '- + ~^ M - <13) 

Linearizing both sides in 9 and Ax and noting that the constant terms in the two linearizations 
match, since they give the on-axis equation satisfied by g'(x), we now get after employing Corol- 
lary 1 that 

,^. (1± ^«.g) ll (14) 
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Moving Ax terms to the right-hand side and simplifying, we get 

~ 1 + g"(x) (S(x,x)+g(x) - h(x)) 
S(x, x) 



b{x, x) 

Differentiating the expression for g'(x) given in Corollary 1, we get that 

dx/dx-l A(x)-1 
9 ( x ) = 5 = 5 • ( 16 ) 

-TO -TO 

Substituting (16) into (15), we get the claimed first-order relationship between 9 and Ax. The 
relationship between 9 and 9 then follows from Theorem 1. □ 



From these expansions, we can compute the phase of the wavefront at the exit pupil or, equiv- 
alently, the path length from an isophase line at the entrance pupil through the system to the exit 
pupil. To this end, note that, for any y , the line y = y + x sin 9, x e [b, c], is an iso-phase line 
for the entering wave. Fix such a y . Let e(x, 9) denote the length of the light ray from this line 
to the first mirror surface. The interval [— f , f ] of the a>axis can be regarded as an exit pupil. Let 
x denote the position along this pupil where the off-axis light beam exits the system. Let f(x, 9) 
denote the length of the light ray from (x + Ax, g(x + Ax)) to this point on the pupil. The path 
length from the iso-phase line to the exit pupil is 

P(x) = e(x, 9) + S(x, x + Ax) + f(x, 9). (17) 



Theorem 2 P(x) =P + y + x9 + o(9). 



Proof. Simple geometry reveals that 

e(x, 9) = y — h(x) cos 9 + x sin 9 = y — h(x) + x9 + o{9) (18) 

and that 

faff) = 9(x + Ax) = g{3t)+g/{3t)M + o{M) 
cos 9 

= g(x)+g'(x)S(x,x)9 + o(9). (19) 



We need a first-order expansion of S(x, x + Ax): 

r,/ ~ a ~x n, ~\ x — x + (q(x) — h(x))q'(x) . „ , . 

S(x, x + Ax) = S(x, x) + — . v " y ; Ax + o(Ax). (20) 

<S(x, x) 

From Corollary 1 and equation (5), we get that 

~ x _ x + (^( 5 ) _ h(x))</(x) = (-P + 5(5) " &(a0) ff'(x) = -S(x, x)</ (x). (21) 
Substituting this into (20) and using Theorem 1, we get 

S(x, x + Ax) = S(x, x) - g\x)S(x, x)6 + o{6). (22) 
Adding terms and using (5), we get the claimed result. □ 



Lemma 3 x = x + (S(x, x) + g(x)A(x)) 6 + 0(6). 

Proof. Simple geometry followed by application of previously derived first-order expansions 
reveal that 

x = x + Ax + g(x + Ax) tan# 

= x + Ax + (g(x) + g'(x)Ax)9 + 0(6) 

= x + (S(x,x)+g(x)A(x))9 + o(6). (23) 
This completes the proof. □ 



Now consider a focusing element placed over the interval [— §, f] of the positive x-axis of 
focal length / as shown in Figure 1 . The electric field at the image plane is 

b 

E(0 = J e' l <^- P(x)+Po+yo hx. (24) 
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where k = 2ir/\ is the wave number. Substituting the first order expansions for x and P(x), we 
get 



b 



E(£) = / e~ lk[ - f ' > M >dx 



d 
2 

/ 

_ d 
2 



■■ / x + (S(x(x),x) + g(x)A(x))0 f r (z.)q\ . 

l H / « xla,;tf M(a;)dx. (25) 



The second integral follows from the substitution x = x(x). 

3.1. Recasting in Unitless Terms 

To write the electric field using unitless quantities, we introduce 6, 6, and x defined as 
follows: 

e = e±, ~e = e- d , £ = £^, x = xd (26) 

where D = b — a denotes the aperture of the primary mirror. Each of these quantities is unitless. 
Associated with these changes of units, we introduce the following normalized versions of the 
apodization function, the transfer function, etc.: 

A(x) = A(xd)d/D (27) 

x(x) = x(xd)/D (28) 

S(x) = S(x(xd),xd)/D (29) 

g(x) = g(xd)/d. (30) 

Making these substitutions, we get 

i 

2 



The term in the exponential with the X/d factor is orders of magnitude smaller than the other terms, 
hence we drop it. The result is an explicit expression for the electric field in the image plane as a 
function of incidence angle 6, which is valid for small 6: 
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1 

2 

Theorem 3 To first order in 0, E(£) = D J e - 2m ^- x ^ A(x)dx. 



It is easy to check that the normalized apodization function is related to the normalized transfer 
function in the same way as before normalization: 



Also, the domain of the A(-) is [— \, \\ and it integrates to one. 



(32) 



The magnitude of the complex electric field is easy to compute if we assume that A(-) is an 
even function. With this assumption, it follows that 



X 

x(x) - £c(0) = J A(u)du 



(33) 



is an odd function. Hence, 



i 

2 

E{£) = De 27rix{0)e J cos (2tt (x£ - (x(x) - x(0))0)) A(x)dx. (34) 



The point-spread function (psf), which is the square of the magnitude of the electric field, is 
then given to first order by 

2 



Psf(£) = D 2 



\ 



cos (2tt (x£ - (x(x) - x(0))0)) A(x)dx 



(35) 



J 



3.2. Nonlinear Effects 



If the expression x(x) — x(0) is linear, then the off-axis psf is just a shift of the on-axis psf. 
Deviation from linearity, on the other hand, produces image degradation. One measure of such 
deviation is the rms phase error relative to A(0): 

1/2 



<} \ 

/ (x(x)-x(0)- A(0)x) 2 dx 



Phase-Error 



V* 



(36) 
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lt is easy to check that the unitless first-order relationship between input angle 8 and output 
angle is 

6 = A(x)0. (37) 

For a uniform apodization, A(x) = 1 (see next section) and hence the unitless input and output an- 
gles match. However, for a nonuniform apodization, the unitless output angle depends on position 
across the pupil. On average, it equals the unitless input angle. But, the physical reality dictates 
using an amplitude weighted average rather than a simple average. In this case, the average unitless 
output angle is larger than the input angle: 

ii i 

2 2 2 

(0) = J Odx = J OA(x)dx = J A{xfdx > (38) 

_i _i _i 

2 2 2 

(this last inequality follows from the Cauchy-Schwarz inequality and the fact that A(x) integrates 
to one). For the types of apodizations arising in high-contrast imaging, the average unitless output 
angle can be more than twice as large as the input angle (see the numerical example in Section 5). 
This extra magnification is one of the attractive features of this approach to apodization. 



4. Simple Examples 

In this section, we consider a few simple examples. 



4.1. Uniform Apodization 

If the apodization function A(-) is a constant, call it A, then x(x) = x + Ax and its inverse 
is x(x) — (x — x )/A. Plugging these into (9), we get that 

h ' (x) = (1 -*>» + >» (39 ) 

-TO 

9'{x) = ^ + ^x (40) 

-TO -M) 

and hence that each mirror is parabolic. In the special case where A = 1, the two mirrors are 
actually planar and the system is a simple periscope. For A^l, the differential equations integrate 
as follows: 

Hx) = c+<^-^ (41) 

»(*) = c+<£ itr < 42 > 
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where c is an arbitrary constant of integration, £ = —x(0)/(A — 1) is the re-coordinate of the 
centerline of the system, H = AP /(2(A — 1)) is the focal length of the "h" mirror, and G = 
P /(2(A - 1)) is the focal length of the "g" mirror. Note that H = AG and hence that the 
magnification of the system is A. Using the fact that A = D/d, the normalized apodization 
function is 



A(x) = 

and the normalized transfer function is 

x(x) 



A(xd) — = Ad jD = 1 



x(xd) xq 



+ X. 



(43) 



(44) 



D D 

Since x(x) depends linearly on x, it follows that the off-axis point-spread function is just a shift 
of the on-axis psf: 



Psf(£) = D 2 



( h 

J cos {2nx {£-0))dx 



\ 



V 



(45) 



/ 



Of course, in practice such a pair of parabolic mirrors will exhibit off-axis errors beyond just a 
shift of the psf but these errors can only be captured by a second-order (or higher) analysis. 



4.2. A Cosine Apodization 

A simple explicit apodization function that approximates the sort of smoothly tapering apodiza- 
tions needed for high-contrast imaging is given by the cosine function: 

/ 2nx\ 

A(x) = a I 1 + cos— —J . (46) 

For this case, the transfer function x(x) is given by 

x(x) = x{0) + ax + a — sin — — (47) 

2ir d 

and the "g" mirror's surface is given by 

, x(0)x + 2-^-x 2 - a-f^ cos(27rx) 

g(x) = c + 2 — ^ '-. (48) 

"o 

The expression for the "h" mirror's surface can be given explicitly using (10) and the inverse 
of x(x) but it is messy so we don't record it here. The amplitude parameter a determines the 
relationship between d and D: 

^ = x ^)- x (Q) = a^ =► D = ad. (49) 
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The normalized apodization is 



A(x) = 1 + cos(27rx) 



(50) 



and the normalized transfer function is 



x(x) = — + x + — sm(2nx). 



(51) 



The nonlinearity of the transfer function implies a nontrivial transformation in addition to the usual 
shift of the off-axis psf: 



In the next section, we consider a similar smooth apodization and plot some of the off-axis 
psf 's to see explicitly the impact of this nonlinear effect. 



The apodization function shown in Figure 3 has unit area and provides 10~ 10 contrast from 
AX/D to QOX/D. It was computed using the methods given in Vanderbei et al. (2003). The on-axis 
point spread function, computed by setting 9 = in (35) is shown in Figure 4. 

In a solar-like system, a star's diameter is about O.Olau. Hence, if a lau Earth-like planet 
appears from Earth at 5X/D, then the angular extent of the star would be from — 0.025A/-D to 
0.025A/-D. Therefore, starlight will enter the system from angle 6 = 0.02X/D. The psf associated 
with this slightly off-axis starlight is shown in Figure 5. Although it is not possible to see it from 
the figure, the threshold for high contrast, i.e. 10 -10 , is pushed to 6.0X/D. Thresholds for high- 
contrast associated with other values of are shown in Table 1. It is clear from this table that 
non-zero stellar size and pointing error increase the inner working angle to at least 6X/D. 

The psf's associated with a planet at 5 and 10 X/D are shown in Figures 6 and 7. Note that 
the magnification due to nonlinearity, described at the end of Section 3, is clearly evident. For 
example, in Figure 6, the input angle is 5X/D but the peak of the psf occurs at about 12X/D. 
Hence, there is a magnification by more than a factor of two. 

On the other hand, the unsharpening of the psf's and the corresponding reduction in peak value 
suggest that a contrast level of 1CT 10 for the on-axis case might not be sufficient for planet detection 
using only the two-mirror system discussed here. Guyon (2003) suggests that an additional set of 




(52) 



5. A Numerical Example 
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two mirrors can be used to restore the off-axis (planet) images with good fidelity. The detailed 
investigation of this addition to the optics will be investigated in a future paper. 
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2.5 




Fig. 1. — Parallel light rays come down from above, reflect off the bottom mirror (on the left), 
bounce up to the top mirror (on the right), and then exit downward as a parallel bundle with a 
concentration of rays in the center of the bundle and thinning out toward the edges — that is, the 
exit bundle is apodized. A focusing element on the rr-axis then creates an image. 
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Fig. 2. — Notations for off-axis expressions. Parallel light rays come down at an angle 6 from 
vertical, reflect off the bottom mirror (on the left), bounce up to the top mirror (on the right), and 
then exit downward at an angle 9 hitting the rr-axis at a new location x. 
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Fig. 3. — A unit area apodization providing contrast of 1CT 10 from A\/D to 60A/D. 




-160 - 



, 1 1 1 1 1 1 

10 20 30 40 50 60 



Fig. 4. — The on-axis point spread function for the apodization shown in Figure 3. 
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Fig. 5. — The off-axis point spread function for the apodization shown in Figure 3 computed at 
9 = OMX/D. 
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Fig. 6. — The off-axis point spread function for the apodization shown in Figure 3 computed at 6 = 
5X/D. Note that the main lobe appears at 12X/D. The nonuniformity of the apodization accounts 
for both the spreading out of the point spread function as well as the apparent magnification by a 
factor of 12/5 = 2.4. 



-20- 




-60 -40 -20 20 40 60 



Fig. 7. — The off-axis point spread function for the apodization shown in Figure 3 computed at 
9 = 10X/D. Note that the main lobe appears at about 22X/D for an effective magnification of 
about 2.2. 
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Table 1 : Output image properties as a function of input angle. Column 1 is the input angle on the 
sky of a point source. Columns 2 and 3 are the output angles at which the intensity in the focal 
plane falls to 1CT 10 compared to the on-axis image peak. Column 4 is the output angle at which 
the star image has its peak intensity. Column 5 is the effective angular magnification factor, i.e., 
the ratio of column 4 to column 1 . 



